Error- and Loss- Tolerances of Surface Codes with General Lattice Structures 
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We propose a family of surface codes with general lattice structures, where the error-tolerances 
against bit and phase errors can be controlled asymmetrically by changing the underlying lattice 
geometries. The surface codes on various lattices are found to be efficient in the sense that their 
threshold values universally approach the quantum Gilbert- Varshamov bound. We find that the 
error-tolerance of surface codes depends on the connectivity of underlying lattices; the error chains 
on a lattice of lower connectivity are easier to correct. On the other hand, the loss-tolerance of surface 
codes exhibits an opposite behavior; the logical information on a lattice of higher connectivity has 
more robustness against qubit loss. As a result, we come upon a fundamental trade-off between 
error- and loss-tolerances in the family of the surface codes with different lattice geometries. 



Introduction. — Recently, topological order has at- 
tracted much interest in both condensed matter physics 
[l|, [H and quantum information science Q . The ground 
state degeneracy of topologically ordered phase cannot 
be distinguished by local operations and hence robust 
against local perturbations. By encoding quantum infor- 
mation into such topologically degenerate subspaces, so- 
called topological quantum error correction (QEC) codes 
B El i logical information can be protected from decoher- 
ence by repeated quantum error correction. There have 
been two types of the topological QEC codes, so-called 
surface codes Q and color codes Q , both of which are the 
CSS (Calderbank-Shor-Steane) codes @. In both cases, 
the threshold values under perfect syndrome measure- 
ments have been calculated to be ~ 11% [1,0], which is 
close to the quantum Gilbert- Varshamov bound [B| in the 
limit of zero asymptotic rate with symmetric X and Z 
errors. In the case of the color codes, their performances 
have been compared among different lattice geometries, 
and their thresholds result in similar values ~ 11% @, El- 
This result is reasonable by considering the fact that the 
color codes are self-dual CSS codes, that is, they are sym- 
metric under a Hadamard transformation. The surface 
codes, on the other hand, are not self-dual CSS codes, 
and hence it is possible to break the symmetry between 
properties of X and Z error-corrections. The surface 
code, however, has been intensively investigated so far 
only on the square lattice, which is a self-dual lattice, 
and therefore its error-correction properties are symmet- 



In this letter, we investigate the surface codes with 
general lattice geometries. Their constructions and error- 
correction procedures are basically the same as those of 
the original surface code. Since the stabilizer operators 
are not always symmetric under the duality transforma- 
tion of the lattice (i.e. exchange of the vertexes and faces 
with each other), the error-tolerances of the surface codes 
with general lattice geometries arc not always symmetric 
between X and Z errors. Interestingly, we find that such 



asymmetry in the error-tolerance is related to the con- 
nectivity of the lattice which defines the surface; error 
chains on a lattice of lower connectivity can be corrected 
easily. Intuitively, this can be understood that finding 
appropriate pairs of incorrect error syndromes, which are 
the boundaries of the error chains, on the lattice of lower 
connectivity is easier, since the incorrect error syndromes 
are more isolated and less percolative. Furthermore, we 
find that the threshold values for independent X and Z 
errors exhibit a universal behavior; they, independent of 
the lattice geometries, approach the quantum Gilbcrt- 
Varshamov bound [H[ in the limit of zero asymptotic rate 
with asymmetric X and Z errors. In this sense, the fam- 
ily of the surface codes can be said to be efficient. We 
also provide a recursive way to construct highly asym- 
metric surface codes on fractal-like lattices. In many ex- 
perimental situations, dephasing is a dominant source of 
errors and therefore the present family of asymmet- 
ric surface codes will help us to correct such biased noise 
efficiently. 

The loss-tolerant scheme is based on a bond percola- 
tion phenomenon, where a reliable logical operator can 
be reconstructed even on the lossy surface as long as the 
survival probability of the qubits is higher than the bond 
percolation threshold [lfj| . Thus the logical information 
on a lattice of higher connectivity is robust against qubit 
loss. As a result, we come upon a fundamental trade-off 
between error- and loss-tolerances of the surface codes 
depending on the connectivity of the underlying lattices; 
the logical information on a lattice of higher connectivity 
is robust against qubit loss, but the error chains on such 
a lattice are difficult to correct (and vice versa). It is in- 
teresting to note that such a trade-off between error- and 
loss-tolerances has been also discussed in a far different 
situation [111 ]. 

Surface codes on general lattices. — Let us consider a 
lattice C(E, V, F), where E, V, and F are the sets of 
edges, vertices, and faces of the lattice, respectively [see 
Fig. Q] (a)]. A qubit is associated with each edge. The 
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FIG. 1: (Color online) (a) A surface code with a general 
lattice structure, (b) The surface code on the hexagonal lat- 
tice and the logical operators, (c) Incorrect error syndromes 
against Z (left) and X (right) errors, which are associated 
with the chains on the primal and dual lattices respectively. 
Here, the Z (left) and X (right) errors are located at the 
same qubits depicted by the solid circles, (d) The remaining 
edges (qubits) on the primal (left) and dual (right) lattices 
after qubit loss. The logical Z and X operators are associ- 
ated with the non-contractible loop on the primal and dual 
lattices, respectively. 



stabilizer operators of the surface code on the lattice £ 
are defined for each face / G F and vertex v £ V, respec- 
tively, as 
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Here, Zi and Xj denote Pauli operators on the ith 
and jth qubits S E) respectively, and Ef and E v 
(Ef )V C E) indicate the sets of edges which are surround- 
ing the face / and are adjacent to the vertex v, respec- 
tively. From these definitions we define the lattice L as 
primal. Its dual lattice, which is denoted by C(E, V, F), 
can be defined by exchanging the vertices and faces of C 
with each other. Similarly to the original surface code 
on the square lattice, Z and X errors are detected as in- 
correct error syndromes at boundaries of Z and X error 
chains on the edges of primal and dual lattices, respec- 
tively. The logical Z and X operators arc defined as 
tensor products of Pauli Z and X operators on the non- 
contractible loops on the edges of the primal and its dual 
lattices [see Fig. [2(b)], respectively. In the following, the 
logical operators play the following two roles: (i) If the er- 
ror correction is failed, a logical operator acts on the code 
space, which we call a logical error, (ii) The logical oper- 
ators represent encoded quantum information, which we 



call logical information. Since H®WAfH®W = B$, the 
surface code defined on the primal lattice C is equivalent 
to that defined on the dual lattice C up to the Hadamard 
transformation iJ®l- E l. In the case of the original surface 
code, which is defined on the square lattice, the Z and 
X error-tolerances are symmetric, since C = C up to 
translation. 

Let us consider, for instance, a surface code on the 
hexagonal lattice as shown in Fig. Q] (b) . In such a case, 
C 7^ £, and hence the Z and X error-tolerances are no 
longer symmetric. Suppose that Z or X errors are, for 
example, located at those qubits depicted with solid cir- 
cles in Fig. [T](c). The Z and X errors are detected at the 
boundaries (vertices) of the error chains on the edges of 
the hexagonal (primal) and triangular (dual) lattices, re- 
spectively, which are denoted by blue (left) triangles and 
red (right) hexagons in Fig. [T](c). In the error-correction 
procedure, we have to infer the locations of the errors 
from the information of incorrect error syndromes. As 
can be seen clearly from Fig. Q] (c), X errors (right) are 
less distinguishable than Z errors (left), since incorrect 
error syndromes are more dense and percolative. More 
precisely, the error syndromes of Z and X errors are asso- 
ciated with the vertices on the primal and dual lattices, 
respectively, and therefore the distinguishability of the 
Z and X errors depends on the connectivity of the pri- 
mal and dual lattices, respectively. The bond-percolation 
threshold is one of the most well-known measures of the 
connectivity of lattices. The bond-percolation thresh- 
olds of the hexagonal (primal) and triangular (dual) lat- 
tices are given by 0.653 and 0.347, respectively. Since 
the hexagonal (primal) lattice has lower connectivity, the 
surface code on the hexagonal lattice is expected to be 
more robust against Z errors. 

Threshold values. — In order to confirm the above ob- 
servation, we perform numerical simulations and esti- 
mate the threshold values of the surface codes with var- 
ious lattice geometries such as the Kagome, hexagonal, 
and (3, 12 2 ) (triangle-hexagonal), say tri-hexa, lattices. 
For independent bit and phase errors with probability px 
and pz, respectively, topological error correction is simu- 
lated by using the minimum-weight-perfect-matching al- 
gorithm. Then, by using the finite-size scaling ansatz 
for the logical error probability, p l z x = A + B(pz,x — 
p t z X )L 1 / 1 ' with fitting parameters (A, B, v, p^x ) s i m ~ 
ilarly to Ref. 12j, the threshold values p^ and p^ 



agains t phys ical Z and X errors arc calculated, respec- 
tively 13j,[l4|. The resultant threshold values are summa- 
rized in Table |U As we expected, when the lattice is not 
self-dual, the threshold values of the surface code exhibit 
asymmetry between Z and X errors, and their behavior 
under various lattices can be understood from the con- 
nectivity of the underlying lattices; if a primal lattice has 
lower connectivity (i.e. higher bond-percolation thresh- 
old), the corresponding surface code is robust against Z 
errors (i.e. higher threshold value p^). On the other 
hand, in such a case, the dual lattice exhibits higher 
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lattice 


Pz 


bond-percolation (primal) 


Px 


bond-percolation (dual) 


square 


0.103 


1/2 


0.103 


1/2 


Kagome 


0.116 


0.524 


0.095 


0.476 


hexagonal 


0.159 


0.653 


0.065 


0.347 


tri-hexa 


0.205 


0.740 


0.041 


0.260 



TABLE I: Summary of the threshold values of the surface codes on various lattices. The bond-percolation thresholds of the 
primal and dual lattices are also shown. 



connectivity due to the Kesten's duality theorem (i.e. 
the sum of bond-percolation thresholds of a primal and 
its dual lattices is equal to one) 15j. Thus the thresh- 
old value of the X errors becomes lower [l6|. Although 
there is no rigorous correspondence between the bond- 
percolation model and the topological error correction, 
it can be used as a good measure of the performance of 
the surface code. 

In Fig. [21 these threshold values (pjf'Pz 1 ) arc plot- 
ted with the quantum Gilbert- Varshamov bound R = 
1 — h(px) — h(pz) in the limit of zero asymptotic rate 
R —> 0, where h{p) is the binary Shannon entropy In- 
terestingly, for all lattices considered here, the threshold 
values {PxiPz) exhibit a universal behavior: they, inde- 
pendent of the lattice geometries, approach the bound. 
In this sense, the present family of the surface codes can 
be said to be efficient. Furthermore, the scaling expo- 
nents are obtained to be v 1.5 for all lattices, which in- 
dicates these surface codes belong to the same universal- 
ity class as that on the square lattice 12[ . From the above 
observation and these numerical evidences, it is expected 
in general that there is a trade-off between X and Z error- 
tolerances of the surface codes with other lattice geome- 
tries, and the trade-off is well characterized by the quan- 
tum Gilbert- Varshamov bound. Since the minimum- 
weight-perfect-matching algorithm is a sub-optimal de- 
coding method, the threshold values will be improved by 



using more sophisticated decoding algorithms |10L Il7| 

There is an exact correspondence between the thresh- 
old values p of the surface code and the critical point 
of the random-bond Ising model on the Nishimori line 
e -2^ c = p th /(l — p th ) where j3 c being the critical inverse 
temperature [3j. Thus, if the threshold value p^ for the 
Z errors is higher, the corresponding random-bond Ising 
model on the primal lattice exhibits a higher critical tem- 
perature. On the other hand, in such a case, the thresh- 
old value Px for the X errors is lower, and hence the 
critical temperature of the random-bond Ising model on 
the dual lattice becomes lower. According to the above 
result, the trade-off between the critical temperatures of 
the primal and dual lattices are also characterized by the 
Gilbert- Varshamov bound through the Nishimori tem- 
perature e~ 2l3c = p th /(l — p th ). 

Loss-tolerance. — Next we consider the performance of 
the surface codes against qubit loss. Recently, a novel 
loss-tolerant scheme has been proposed for the surface 




FIG. 2: (Color online) The threshold values (px,Pz) for the 
square, Kagome, hexagonal, and tri-hexa lattices are plotted 
(red dots) with the quantum Gilbert- Varshamov bound (green 
line). The blue dots are the threshold values for the surface 
codes defined on the duals of the Kagome, hexagonal, and 
tri-hexa lattices (from left to right). 



code on the square lattice [lfjj. There, even if some qubits 
are lost, the stabilizer operators and the logical operators 
can be reconstructed as long as the loss rate is less than a 
threshold. The threshold value of qubit loss without er- 
rors is given exactly by the bond-percolation threshold of 
the lattice. This is because one can reconstruct a logical 
operator as long as the remaining qubits (bond) percolate 
through the lossy surface [see Fig. Q] (d)]. We can also 
apply the loss-tolerant scheme for the surface codes on 
the general lattices in a similar manner as in Ref. fioj ]. 
Then, the loss-tolerances of the logical X and Z infor- 
mation also become asymmetric depending on underling 
lattice geometries. According to the Kesten's duality the- 
orem, there is a trade-off between the bond-percolation 
thresholds of the primal and dual lattice. Thus, asym- 
metry in the tolerable loss rates p x ss and p^ ss of the 
logical X and Z information, respectively,is subject to a 
trade-off, p x ss +p'| ss = 1. 

We also perform numerical simulations for the loss- 
tolerant scheme on the Kagome, hexagonal, and tri-hexa 
lattices. The resultant threshold curves (p loss ,p th ) of the 
error and loss rates are plotted in Fig. [3] (a) for each 
lattice for logical X and Z information. The logical X 
information is robust against both Z errors and qubit loss 
as shown in Fig. [3] (a) , while the logical Z information is 
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fragile against both X errors and qubit loss |l6[. These 
results can be understood from the two trade-offs, which 
we have studied above, the trade-off between the Z and 
X error-tolerances [points on the vertical axis in Fig. [3] 
(a)], and the trade-off between the loss-tolerances of log- 
ical Z and X information [points on the horizontal axis 
in Fig. [3] (a)]. If we choose the primal lattice C of low 
connectivity so that the surface code is robust against Z 
errors [see Fig. Q] (c) left] , then the bonds of such a pri- 
mal lattice are less percolative, and hence the logical Z 
information becomes fragile against loss [sec [T] (d) left] . 
This is also the case for fragility against X errors and 
robustness of the logical X information [sec Fig. [T] (c) 
right and (d) right]. That is, there is a fundamental rela- 
tionship between robustness against errors and fragility 
of the logical information against qubit loss depending 
on the lattice geometries of the surface codes. 

Systematic construction of asymmetric lattices. — In 
the above case, qubit loss occurs randomly according to 
the loss rate p loss . It is also suggestive to consider a sys- 
tematic qubit loss. If we introduce a systematic qubit loss 
for the triangular (dual of the hexagonal) lattice, whose 
threshold value is 15.9% as shown in Table HI then we 
can reconstruct the surface as the square lattice, whose 
threshold value is 10.3% [sec Fig. [3](b) C and &]. By 
doing the inverse operation of the systematic qubit loss, 
namely, a systematic injection of qubits, we can construct 
a new lattice C n recursively, which has high connectivity 
(its dual lattice C n has low connectivity), as shown in 
Fig. [3] (b). Therefore, the threshold value for the Z (X) 
errors on C n (C n ), which has lower (higher) connectiv- 
ity, is expected to be higher (lower) than that of £„_i 
(£„_i). On the other hand, the logical Z (X) informa- 
tion is more fragile (robust) against qubit loss. 

Conclusion and discussion. — We have investigated the 
surface codes with general lattice geometries. These sur- 
face codes are efficient, since their threshold values ap- 
proach the quantum Gilbert- Varshamov bound univer- 
sally. We have found that there are interesting relation- 
ships between error or loss-tolerances of the surface codes 
and the geometrical properties of the underlying lattices, 
where the connectivity of the lattice has an opposite ef- 
fect in error and loss-tolerances. This results in a funda- 
mental trade-off between robustness against errors and 
fragility of the logical information against qubit loss (and 
vice versa) among the surface codes on the different lat- 
tices. We also provide a recursive way to construct a 
highly asymmetric surface code on a fractal-like lattice. 
These results suggest that we can systematically find an 
efficient surface code by selecting the underlying lattice 
geometry depending on the degree of biased noise. 

The surface codes with general lattice geometries 
also support topologically protected measurement-based 
quantum computation by braiding the defects (Til l20j ]. 
The error- and loss-tolerances in fault-tolerant topologi- 
cal quantum computation on the general lattices are an 
intriguing future problem. 
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FIG. 3: (Color online) (a) The threshold curves (p loss ,p th ) for 
logical X (solid line) and Z (dotted line) information [Tg] ■ Top 
to bottom: tri-hexa (X), hexagonal (X), Kagome {X), square 
(X and Z), Kagome (Z), hexagonal (Z), and tri-hexa (Z). (b) 
A recursive way to construct a fractal-like lattice C n of higher 
connectivity. Its dual lattice £ n has lower connectivity. 



It is an interesting correspondence that a trade-off be- 
tween error- and loss-tolerances has been found also in a 
far different situation [ll[, where lost information is re- 
covered by using counter-factual indirect measurements 
on a specific shape of the cluster states. In such a case, 
if we introduce more redundancy for the loss-tolerance, 
more errors are accumulated, which results in a similar 
trade-off. This correspondence suggests that there might 
be more fundamental principle between error- and loss- 
tolerances; for two mutually noncommuting observables, 
if one observable is robust against qubit loss, the other 
one must be fragile against environmental perturbations 
21 1 . 

KF is supported by MEXT Grant-in- Aid for Scientific 
Research on Innovative Areas 20104003. 

Note added. — During preparation of this manuscript 
[22} , a related work has been appeared 
supports generality of the present results. 
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